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We report preliminary results of our study of matrix elements of baryon number violating operators which 
appear in the low-energy effective Lagrangian of (SUSY-)Grand Unified Theories. The calculation is performed 



on a 32'' x 80 lattice at f3 



.1 using Wilson fermions in the quenched approximation. Our calculation is 



independent of details of (SUSY-) GUT models and covers all interesting decay modes. 



1. Introduction 

Proton decay is one of the most exciting predic- 
tions of Grand Unified Theories (GUTs). Experi- 
mental effort over the years has pushed the lower 
limit on the partial lifetimes to t / Bp^^o^+ > 
5.5 X 10^2 years and r/Bp^j^+i; > 1.0 x 10^2 
years at the 90% confidence level||l[], and fur- 
ther improvement is expected from the Su- 
perKamiokande experiment. 

A crucial link to relate these numbers to con- 
straints on (SUSY-)GUT models is the values of 
hadron matrix elements relevant for proton decay. 
Model calculations suffer from high degree of un- 
certainty, various estimations easily differing by a 
factor ten. 

Pioneering lattice QCD studies to remove this 
source of uncertainty were made about ten years 
ago, first combining calculations of the matrix 
element {0\O^\p) and soft-pion theorems to es- 
timate (tt" I |p) |^,|| , and subsequently directly 
evaluating the matrix element {■k'^\0^\p) itself 

In this article we report preliminary results of 
our renewed effort to determine the matrix ele- 
ments from first principles of QCD. In addition 



to the use of larger lattice sizes and higher statis- 
tics to achieve a much better precision, which is 
made possible through the increase of computing 
power, we aim to advance the calculation on two 
fronts: (i) calculation of p K matrix elements 
relevant for SUSY-GUT as well as those for the 
p TT mode, for physical values of u-d and s 
quark masses and for physical momenta, and (ii) 
evaluation of matrix elements of all dimcnsion- 
6 baryon number violating operators classified 
according to SU{3)c x SU{2)l x U{1)y invari- 
ance[^|j, so as to cover various GUT models and 
decay processes. 



2. Calculational procedure 

Our calculation is carried out in quenched QCD 
at (3=6.1 with the Wilson quark action on a 
32^ X 80 lattice. We analyze 100 configurations 
for the hopping parameter i^=0. 15428, 0.15381, 
0.15333, 0.15287. The lattice scale fixed by rup = 
769MeV in the chiral hmit {K^ = 0.15499(2)) 



equals 



'presented by N. Tsutsui 



2.56(4)GeV, and the point for 
= 0.648 is 

given hy K = 0.15304(5). All errors are estimated 
by the single elimination jackknife procedure. 
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Figure 1. Matrix element relevant for p ^ tt'^ 
and p — > K'^ decay normalized by measured pro- 
ton mass at p = of pseudo scalar meson as a 
function of u-d and s quark mass in lattice units. 
u-d quark mass is taken to the chiral limit for Ak- 

To calculate the nucleon (N) to pseudo scalar 
(PS) meson matrix element of a baryon number 
violating operator , we form the ratio, 
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{PS\0^\N) = 



{0\JpsO^ Jn\0) 
{0\JnJn\0){0\Jps4s\0) 



X {0\Jps\PS){N\Jn\0), 



(1) 



where {0\Jps\PS) and {N\Jn\0) are extracted 
from local-local hadron propagators. We fix the 
nucleon source at t=0, PS sink at t = 32 and move 
the operator between them. Matrix elements are 
evaluated for four spatial momenta pa=(0,G,0), 
(1,0,0), (0,1,0), (0,0,1) in units of minimum mo- 
mentum ap„iin = 7r/16 injected in the PS meson. 

We distinguish u-d and s quark masses; the 
former is taken to the chiral limit, and the latter 
interpolated to the physical s quark mass in our 
calculations. After this procedure, we interpolate 
the spatial momentum to the physical value. 

Matrix elements are renormalized, with mixing 
included, by tadpole-improved one-loop renor- 
malization factors to the MS scheme calculated 
at the scale ji = 1/a. 
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Figure 2. Momentum dependence of p ^ 7r° and 
p —)■ matrix elements. Quark masses are taken 
to the physical point. 

3. Results 

Let us define the p — > tt'^ and p — > matrix 
elements A.^ and Ak by 

(^°(/^°)|e.,,(^^C4(4))4b) EE A.(^)iV^. (2) 

In Fig. |l| we show our results for these matrix 
elements at zero momentum p = 0, normalized 
by proton mass measured for relevant values of 
u-d quark mass. For A^r, abscissa represents the 
u-d quark mass, while it represents the s quark 
mass for Ak, the u-d quark mass having been 
taken to the chiral limit. 

Compared to a first calculation of the Atj ma- 
trix element [Q, whose results are consistent with 
ours, our statistical errors of 4-8% are improved 
by about a factor five. This allows us to ob- 
serve that the amplitude exhibits a clear de- 
crease of about 40% from the region of s quark 
(jUu^dO' ~ 0.04) to the chiral limit, which was not 
apparent in results of Ref. Q| . 

For the p K'^ matrix elements Ak, after chi- 
ral extrapolation of u-d quark mass, the depen- 
dence on the s quark mass is small. The point 
plotted with a cross shows the value interpolated 
to the physical s quark mass. 

In Fig. ^ we plot A^^ and Ak as a function of 
squared momentum after quark masses are taken 
to the physical values. We observe a quite signifi- 
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Figure 3. Matrix elements for p n'^ , tt'^ , K'^ 
decay modes from present work (open circles) 
compared with predictions of tree-level chiral La- 
grangian (asterisks). 

cant momentum dependence for both amplitudes, 
necessitating an interpolation for a precise esti- 
mate of their physical values. The open circles in 
Fig. H show results of a linear interpolation in 
to physical momentum. 

Phenomenological analyses of proton decay of- 
ten employ the soft-pion relation 



V2U 



(l + F + D) 



(3) 



to estimate the p 



matrix element, where F 



and D arc the axial vector matrix elements of pro- 
ton. With our results for physical quark masses 
and momentum, the right hand-side is about four 
times larger than the left hand-side of this re- 
lation. A similar discrepancy was observed in 
Ref. 01 . Examining the above relation for zero 
pion momentum in the chiral limit, i.e., the real 
soft-pion limit, we find that the two sides are still 
discrepant by about a factor two. The origin of 
the discrepancy is not clear to us at present. 

4. Phenomenology 

We plot our preliminary results for matrix el- 
ements normalized by proton mass squared rel- 
evant for p 7r°,7r+ and K'^ decay in Fig.|[ 
For comparison, predictions from tree-level chi- 



ral Lagrangian with a choice of the parameters 
—a ^ 13 ~ 0.003 GeV'^ are shown, which rep- 
resent the smallest values among various model 
estimations. Here a and (3 are defined as 
(0|eyfc(w'C4(^))M^|p) = a{l3)NR. Our lattice 
results are even smaller for some of the chan- 
nels, and predictions for the ratio of p — ^ tt" and 
p amplitudes are also different. 

For non-SUSY minimal SU(5) GUT, the decay 
width for the p ^ ti^ mode is given 



T{p 



) 



\-al{M^)Al^, 



(4) 



X 



Employing the GUT gauge coupling a5[Mx) = 
0.024 and the short distance renormalization fac- 
tor = 10 at ^ = 2GeV as in Ref. [|, 
substituting our preliminary result for At^ yields 



^p^7r"e+ V-^--^ -1- u.-l; ^ ^kj V2.0xlO"GeV 

years where the error is only statistical. 



= (1.1 ± 0.1) X 10^"( 
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